CONSTRUCTION OF THE PATHS OF BROWNIAN MOTIONS 
ON STAR GRAPHS 

VADIM KOSTRYKIN, JURGEN POTTHOFF, AND ROBERT SCHRADER 

Abstract. Pathwise constructions of Brownian motions wliich satisfy all possible 
boundary conditions at the vertex of star graphs are given. 



1. Introduction and Preliminaries 

In the recent years there was a growing interest in metric graphs because of their 
wide range of important applications, see, e.g., the articles in fl^ and the references 
given there. The simplest metric graphs are star or single vertex graphs: They can be 
defined as a set having a finite collection of subsets isomorphic to M+, called external 
edges, where the points corresponding to the origin of ]R_|_ under these isomorphims 
are identified and form the vertex of the graph. One may visualize a star graph as a 
finite number of rays in going out from the origin. 



On the other hand, in his pioneering articles 1 11 - 13 1 Feller investigated the prob 



lem of characterization and construction of all Brownian motions on intervals. This 



problem found a complete solution in the work of Ito and McKean 1 19 20|. In par- 
ticular, in | |19J Ito and McKean solved the problem of construction of the paths of 
all Brownian motions on the semi-line M+ employing the theory of the local time of 



Brownian motion pOl, and the theory of (strong) Markov processes ||3]|6]-[8 16 1. 



Therefore it is natural to investigate Feller's problem on metric graphs, and in par- 
ticular on star graphs. In fact, the Walsh process introduced by Walsh in | ^38J as 



a generalization of the skew Brownian motion |20| is the most basic example of a 
Brownian motion on a star graph, and in the present article it will serve — together 
with its local time at the vertex — as the main building block of our constructions. 
The Brownian motions constructed here on star graphs are then the basic building 



blocks of Brownian motions on general, finite metric graphs in the article 1 24 1 of the 
present authors. 

The main ideas for the construction of Brownian motions with boundary condi- 
tions at the vertex compatible with Feller's theorem (see below, theorem [T3] ) are those 



which can be found in the above mentioned work by Ito and McKean (cf . also 1 23 
Chapter 6]): The reflecting Brownian motion in the case of M_|_ is replaced by a Walsh 
process [38 1 (cf. also, e.g., |j2j) on the single vertex graph, and then slowing down 
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and the killing of this process on the scale of its local time at the vertex are used to 
construct processes implementing the various forms of the Wentzell boundary condi- 
tion. We provide a number of arguments which — at least on a technical level — are 
rather different from those found in the standard literature. For example, whenever 
possible, we use arguments based on Dynkin's formula to derive the domain of the 
generator (i.e., the boundary conditions at the vertex). This approach appears to be 
much simpler and more intuitive than the one with standard arguments |T9|20 23 1 



for the semi-line M+, which is based on rather tricky calculation of heat kernels with 
the help of Levy's theorem. Moreover, for the case of killing, instead of using the 
standard first passage time formula for the hitting time of the vertex we use a first 
passage time formula for the lifetime of the process. In the opinion of the authors this 
leads to much simpler computations of the transition kernels than those in f T9l|20[|23| 
for a Brownian motion on M_|_. In all cases we also derive explicit expressions of the 
analogues of the quantum mechanical scattering matrix on star graphs. 

The article is organized as follows. In several subsections of the present section we 
set up our notation and discuss some preparatory results. In section [2] we recall the 
construction of a Walsh process on a metric graph. In section [3] we construct a Walsh 
process on the single vertex graph with an elastic boundary condition at the vertex, 
while in section [4] we construct a Walsh process with a sticky boundary condition at 
the vertex. The most general Brownian motion on a star graph is obtained in section[5] 

1.1. Brownian on a Star Graph and Feller's Theorem. From now we shall con- 
sider a fixed star graph Q with vertex v and n G N external edges li, . .., In. Q is 
equipped with the natural metric d which is induced by the metric that each exter- 
nal edge inherits from being isomorphic to M+. Thus {Q,d) is a locally compact, 
complete metric space, and we shall always consider Q as equipped with its Borel 
cr-algebra. Q° denotes the set Q \ {v} which we also call — by abuse of language 
— the open interior of Q. It is the disjoint union of n copies = 1, . . . , n, of 
the interval (0, +oo). Every ^ G ^° is in one-to-one correspondence with its local 
coordinates {k, x), where k G {1, . . . , n} is the index of the edge ^ belongs to, and 
X > denotes the distance of ^ to v. For simplicity we shall often write ^ = (/c, x). 

The following definition of a Brownian motion on Q is the generalization of the 
definition of a Brownian motion on the half line R+ as given by Knight in |23| Chap- 
ter 6]. 

Definition 1.1. A Brownian motion X = {Xt, t G M+) on Q is a diffusion process 
on Q, such that X with absorption at v is equivalent to a Brownian motion on the half 
line M-i- with absorption at the origin. 

Remarks 1.2. By a diffusion process we mean a strong Markov process (e.g., in the 
sense of Q), which a.s. has cadlag paths and a.s. the paths are continuous on [0, (), 
where is its lifetime. Moreover, in definition |1.1| we have — as we shall usually 
do without any danger of confusion — identified the external leg 1^, k = 1, 
n, on which the process starts with the corresponding copy of R+. Throughout we 
shall assume without loss of generality that the filtration of a Brownian motion on Q 
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satisfies tlie "usual conditions", i.e., it is rigiit continuous and complete relative to the 
underlying family (P^, G ^) of probability measures. 

Co{Q) denotes the Banach space of continuous functions on Q which vanish at 
infinity equipped with the sup-norm. 

Define Cq{Q) as the subspace of Co{Q) consisting of those functions / G Co{Q) 
which are twice continuously differentiable on Q°, such that /" extends from Q° to 
C/ as a function in Co{Q). The following lemma states some of the properties of 
functions in Cq{Q). It can be proved with the help of applications of the fundamental 
theorem of calculus and the mean value theorem, and the proof is omitted here. 

Lemma 1.3. Suppose that f £ Cq{Q), k G {1, . . . , n}. Then the limit of f'{^) as 
^ converges to v along the open edge exists. The directional derivatives f^^\vk), 
i = 1,2, of first and second order of f at the vertex v in direction of the edge Ik exist, 
and the equalities 

/»(t;fc)= lim i = l,2, 

hold true. Moreover, f vanishes at infinity. 

Remark 1.4. By definition we have that for every / G Cq (^) and all j. A; = 1, . . . , n, 
f"(vj) = f"{vk), and henceforth we shall simply write f"{v) for this quantity. On 
the other hand, in general f'{vj) 7^ f'{vk) for j / k. 

It is not hard to show that every Brownian motion on a star graph is a Feller process. 
A convenient way to prove this is to show that its resolvent maps Co{Q) into itself by 
arguments similar to those in ||20[ Section 3.6], and to observe that the path properties 
imply for all ^ G ^, / G Co{Q), PtfiO converges to f{^) as t decreases to 0, where 
{Pt, t G M+) denotes the semigroup generated by the Brownian motion. Then one 
can use well-known arguments (for example, a complete proof can be found in [27]) 
to conclude the Feller property in its usual form, e.g., [35 1. 



The analogue of Feller's theorem 1 23 Theorem 6.2] for a Brownian motion on the 



single vertex graph Q reads as follows: 

Theorem 1.5. Assume that X is a Brownian motion on Q. Then there exist constants 
a, bk, c G [0,1], k = 1, . . . , n, with 

n 

(1.1a) a + c + ^5fc = l, a /I, 

k=l 

such that the domain T>{A) of the generator A ofX in Co{Q) consists exactly of those 
f G Cl{g) for which 

n 

(1.1b) af{v) + p"{v) = Y,bknvk) 

k=l 

holds true. Moreover, for f e V{A), Af = 1/2/". 
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The proof in (23] for the case where Q has only one external edge, i.e., Q 



is readily modified for a general star graph Q. On the other hand, theorem 1.5 also 



follows from Feller's theorem in the case of a general metric graph [24, Theorem 1.3]. 

1.2. Standard Brownian Motion on the Real Line. The construction of Brow- 
nian motions on a single vertex graph with infinitesimal generator whose domain 
consists of functions / which satisfy the boundary conditions ( |1.1| ) is quite similar 
to the construction carried out for the half-line in f\3l , | [20] |, | [23] |. This in turn is 
based on the properties of a standard Brownian motion on the real line, cf., e.g.. 



|T4J[T5}[T7J|22)|35J|40]|, and the works cited above. For the convenience of the reader, 
and for later reference, we collect the pertinent notions, tools and results here. 

Let {Qx, X G M) denote a family of probability measures on a measurable space 
{Q' , A'), and let B = {Bt, t G M+) denote a standard Brownian motion defined on 
{Q' , A') with Qx{Bq = x) = 1, X G M. It will be convenient to assume throughout 
that B exclusively has continuous paths. Whenever it is notationally convenient, we 
shall also write B{t) for Bt, t > 0. Furthermore, we may suppose that there is a shift 
operator 9 : M+ x O — ^ il, such that for all s,t > 0, Bg o 9t = Bg+t- 

We shall always understand the Brownian family {B, {Qx, x G M)) to be equipped 
with a filtration J' = {J't, t > 0) which is right continuous and complete for the 
family {Qx, x G M). (For example, J' could be chosen as the usual augmentation of 



the natural filtration of B (e.g., \21. Sect. 2.7] or ||35j Sect.'s 1.4, III.2]).) 
For any A C M, we denote by the hitting time of A by B, 

(1.2) = mi{t >0, BtG A}, 

and we note that for all A belonging to the Borel cr-algebra S(M) of M, is a 
stopping time with respect to J' (e.g., [35, Theorem III.2.17]). In the case where 
A = {x}, X G M, we also simply write for H^^y We shall also denote these 

stopping times by H^{A) and H^{x), respectively, whenever it is typographically 
more convenient. The following particular cases deserve special attention. Let x G M. 
Then we have (e.g., (20) Sect. 1.7], 1(22] Sect. 2.6], (35} Sect.'s II.3, III.3]) 



(1.3) Qo{Hx G dt) = Qx{H^ G (it) = ^ g{t,x) dt, t > 0, 
where g is the GauB-kemel 

(1.4) g{t, x) = —L= e-^'/2t^ t > 0, X G M, 

V 2vrt 

and 

(1.5) i?o^(e-"^") = E^ie-^^'o) = g-v^l-l, A > 0. 

Moreover, for a < x < 6 the law of H^^^^ under Qx is well-known (e.g., |20 
Problem 6, Sect. 1.7]), and its expectation is given by 

(1.6) E^{Hf^,^) = {x-a){h-x). 
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Denote by = {Lf , t > 0) the local time of B at zero, where we choose the 
normalization as in, e.g., | |35[ (and which thus differs by a factor 2 from the one used 
in, e.g., | )17^|22J ): for x € M, P^-a.s. 

(1.7) Lf = lim ^A({s<t, <e}), t > 0, 

and here A denotes the Lebesgue measure. Thus, in terms of its a-potential (cf. ||4j 
Theorem V.3.13]) we have 

' e-"*dLf ) = -=e-^l^l, a>0, xGM, 
^ v2a 

which provides an efficient way to compare the various normalizations of the local 

time used in the literature. Slightly informally we can write 

rt 

B 



(1.9) Lf = / 5o{B,)ds, 

Jo 

where 5q is the Dirac distribution concentrated at 0. is adapted to J, and non- 
decreasing. Moreover, for every x S M, Pa;-a.s. the paths of are continuous, and 
is additive in the sense that 

(1.10) Lf+, = Lf + Lf o0i, s,tGM+. 

Similarly as above, we shall occasionally take the notational freedom to rewrite Lf 

as L^{t). 

We will need the following well-known result (e.g., [201 Section 2.2, Problem 3]): 



Lemma 1.6. The joint law of\Bt \ and Lt, t > 0, under Qq is given by 

(1.11) Q^[\Bt\edx,Lf edy)=2^^e-^^^y^"'^'dxdy, x,y>0. 

Let = {Kf, r > 0) denote the right continuous pseudo-inverse of L, 

(1.12) = M{t>0, Lf > r], r > 0. 

Note that due to the a.s. continuity of L^ we have a.s. L^^ = r. In appendix B of 1 25 1 
the present authors proved the following 

Lemma 1.7. For any r > 

(1.13) Qo{Kf e dt) = ^g{t,r)dt, t > 0, 
and 

(1.14) ^0^(6-^^''') =e-^^ A>0 
holds. 

Moreover, we shall make use of the following lemma, which is similar to results in 
Section 6.4 of [22], and which is proved in appendix B of |25 1, too. 

Lemma 1.8. Under Qo, L^ (^H^_^ x > 0, is exponentially distributed with 
mean x. 
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1.3. First Passage Time Formula for Single Vertex Graphs. In this subsection we 
set up some additional notation which will be used throughout this article. Also we 
record a special form of the well-known first passage time formula, e.g., pm[34| . 

Let X be a Brownian motion on Q in the sense of definition | 1 . 1 [ defined on a family 
{'D,,A,T = {Ft, t > 0),{P^, C G G)) of filtered probability spaces. Let be 
the hitting time of the vertex v. It follows from definition |1.1| that for all ^ ^ Q, 
Pj:{Hy < +oo) = 1. For A > 0, set 

(1.15) exiO = E^{expi-XH,)) = e-^^^^'^), ^ e G, 

where E^{-) denotes expectation with respect to P^. The last equality follows from 
formula ( |1.5| ). 

Recall that we denote the natural metric on Q by d. We introduce another symmet- 
ric map dy from Q x Q to M_|_ defined by 

(1.16) dy{C,i]) = d{C,v)+d{v,r,), ^,V^G, 

which is the "distance from to 77 via the vertex v". Observe that if ^, r] £ Q do not 
belong to the same edge, then t]) = d{$^, rj) holds. 
Next we define two heat kernels on Q by 

n 

k=l 
n 

(1.18) py{t, ^,v) = Y. (^)9{t, dy{C, 7?)) h^ir,), 

k=l 

with t > 0, ^, ?7 G t/. 5 is the GauB-kemel defined in equation \\A\ . Hence, in local 
coordinates ^ = (A;, x), rj = (m, y), x, y > 0, k, m e {1,2, . . . , n}, these kernels 
read 

(1.19) p{t,{k,x),{m,y)) = ^e-^^-y^'/^'6km 

V ZTTl 

(1.20) py{t,{k,x),{m,y)) = ^e-^-^^^^' I'-'b^^. 

The Dirichlet heat kernel on Q is then given by 

(1.21) P^{t,C,v)=P{t,C,v)-Pv{t,C,v), t>0,tr]£g. 

It is the transition density of a strong Markov process with state space 5° U { A} which 
on every edge of Q° is equivalent to a Brownian motion until the moment of reaching 
the vertex when it is killed, and A denotes a universal cemetery state for all stochastic 
processes considered adjoined to G as an isolated point. (Observe that this process is 
not a Brownian motion on Q in the sense of definition |1.1[ ) 

The Dirichlet resolvent = [R^ , A > 0) on ^ is defined by 

(1.22) R^f{i) = E^(^j^\-^'f{Xt)dt), x>o,(eg, f eB{g). 
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It is easy to see that has the following integral kernel on Q 

(1.23) rf ??) = rA(e, r?) - r,,A(C, S g, 

where for ^, 77 G ^, 



(1.24) r,(C,r?) = ^ 1^,(0 7=^ hM^ 

k=i 

and 

k=l 

(1.25) 



1 

fc=l 

In particular, is the Laplace transform of the Dirichlet heat kernel ( |1.21[ ) at A > 0. 

In the present context the well-known first passage time formula, e.g., |20j34|, 
reads as follows 

RxfiO = E^(^j\-^'f{Xt)dt)+E^{e-^'^ Rxf{Xs)), 

where S is any P^-a.s. finite stopping time relative to J . The choice S = H^, gives 
the following result. 

Lemma 1.9. Let X be a Brownian motion on Q with resolvent R = {R\, A > 0). 
Then for all X > 0, ^ e G, f G B{g), 

(1-26) Rxf{i) = R^fiO + exiO Rxfiv) 

holds true. 

The following notation will be convenient. For real valued measurable functions 
/, g on Q, with restrictions fk, Qk, k G {1,2, . . . , n}, to the edges 1^ we set 



n 

if, 9) = J f{09{0d^ = Y.^fk,9k), 



k=l 

where the integration is with respect to the Lebesgue measure on Q, and 

{fk,9k)= fk{x)gk{x)dx, 
Jo 

whenever the integrals exist. 

Assume that / G Co{Q). Then for A > 0, R\f belongs to the domain of the 
generator of X, and therefore to Cq{Q) (cf. subsection 1.1 1. It is straightforward to 
compute the derivative of the right hand side of formula ( 1.26l, and we obtain the 
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Corollary 1.10. For every Brownian motion X onQ with resolvent R = {R\, A > 0), 
and all f £ Co{Q), 



(1.27) {Rxfy{vk) = 2{ex,k,fk)-V2\Rxf{v), k e {1,2, . . . ,n}, 
holds true. 

1.4. The Case 6 = 0. The case, where all parameters b^, k = 1, . . . ,n, m equa- 
tion \\.\) vanish, is trivial in the sense that the associated Brownian motion can be 
constructed by a stochastic process living only on the edge where it started, and there- 



fore it is just a classical Brownian motion on M_|_ in the sense of |23 Section 6.1]. 
This case is also discussed briefly in [23| , but for the sake of completeness we include 
it here in somewhat more detail than in fTS\. 

Consider a standard Brownian motion on M as before, and without loss of general- 
ity assume in addition that the underlying sample space is large enough such that all 
constant paths in M can be realized as paths of the Brownian motion. Construct from 
the Brownian motion a new process by stopping it when it reaches the origin of M, and 
then kill it after an exponential holding time (independent of the Brownian motion) 
with rate /3 > 0. We shall only consider starting points x G M+. If /3 = 0, then the 
process is simply a Brownian motion with absorption at the origin. For example, it 
follows from Theorem 10.1 and Theorem 10.2 in [7J that for every /3 > this process 
is a strong Markov process, and obviously it has the path properties which make it a 
Brownian motion on M+ in the sense of fT^, Section 6.1]. Thus, if ^ G t/, ^ G Ik, 
k = 1, . . . , n, then we just have to map this process with the isomorphisms between 
the edges Ik, k = 1, . . . ,n, and the interval [0, +oo) into Q to obtain a Brownian 
motion on Q with start in ^, such that it is stopped when reaching the vertex, and then 
is killed there after an exponential holding time with rate /3 > 0. 

Let = {U^, t > 0) denote the semigroup associated with this process. It is 
obvious that for / G Co{g) we get Uj^f{v) = exp(-/3t)/(t>), t > 0. Thus for the 
corresponding resolvent R^ = (i?^, A > 0), and / G Co(^) one finds 

(1.28) XRxfiv) - f{v) + mlfiv) = 0, A > 0. 



Let be the generator of this process, and recall from theorem 1.5 that for all 
/ G V{A^), ^0/(0 = 1/2 /"(d. ^ e ^- But then the identity Ai?5 = ^"-^a + id 
implies the following formula 

l{Rlf)"{v)+^Rif{v) = 0. 

For every X > R^ maps Co{g) onto With the choice a = (1 + /3)"^ 13, 

c = (1 + this shows that the process realizes the boundary conditions of equa- 
tions ( |1.1[ ) with bk = 0, k = 1, . . . ,n. 



Moreover, we can now use equation ( 1.26 1 combined with formula ( 1.28 1, to obtain 



the following explicit expression for the resolvent with / G Co(^), A > 0: 
(1.29) RlfiO = RxfiO + ^ e-^'^^-'') f{v), e G G, 
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where, as before, is the Dirichlet resolvent. 

In order to compute the heat kernel associated with this process on ^, we invert the 



Laplace transforms in equation ( |1.29[ ). For the first term on the right hand side this is 
trivial, and gives the Dirichlet heat kernel , cf. equation ( |1.21[ ). The second term 
could be handled by a formula which can be found in the tables (e.g., |9, eq. (5.6.10)]). 
But this formula involves the complementary error function erfc at complex argu- 
ments, and does not yield a very intuitive expression. Instead, we can simply use the 
observation that t i— )■ exp(— is the inverse Laplace transform of A i— (/3 + A)^^. 
Moreover, the well-known formula for the density of the hitting time of the origin by 
a Brownian motion on the real line (e.g., 1 ,20 , p. 25], |22, p. 96], [35. p. 102]) provides 
us with the following expression for the density of the first hitting time of the vertex 

(1.30) P^{H, G ds) = e-'^(«''')'/2. , > 0. 
Using the well-known Laplace transform (e.g., |[9j eq. (4.5.28)]) 

(1.31) / e-^' -=e-°-^/'^' ds = e-^'', a > 0, A > 0, 



we infer that the inverse Laplace transform of the exponential in ( |1.29| ) is given by 
P^{Hv G ds). Thus we obtain the following heat kernel 

(1.32) 

p\t,i,dr,)=p^{t,^,7^)d7i-(^j^ e-^(*-^)p^(i/„ G ds)) eMv), 

= p^{t, r,)dv-( f e-/^(*-) ^ e-^(5'^)^/2^ ds) e.{dr,), 
with ^, G ^, t > 0, and is the Dirac measure at the vertex v. 



1.5. Killing via the Local Time at the Vertex. We recall from remark 1.2 that we 
may and will consider every Brownian motion X on ^ with respect to a filtration 
T = {Ft, t > 0) which is right continuous and complete relative to (P^, ^ G Q), and 
such that X is strongly Markovian with respect to F. 

In this subsection we suppose that X is a Brownian motion on the single vertex 
graph Q with infinite lifetime, and such that the vertex is not absorbing. This entails 
(e.g., [35, Proposition II.2. 19]) that X leaves the vertex immediately and begins a 
standard Brownian excursion into one of the edges. Therefore we get in this case 
for the hitting time of the vertex Pv{Hy = 0) = 1, i.e., v is regular for {v} in 
the sense of [4|. Consequently X has a local time L = {Lt, t > 0) at the vertex 
(e.g., [0 Theorem V.3.13]). Without loss of generality, we suppose throughout this 
subsection that L is a perfect continuous homogeneous additive functional (PCHAF) 
of X in the sense of [40, Section III.32]. That is, L is a non-decreasing process, which 
is adapted to F, and such that it is a.s. continuous, additive, i.e., Lt+s = Lt + LgO 6t, 
and for all ^ £ Q, P(,{Lo = O) = 1 holds true. Moreover we may and will assume 
from now on that X and L are pathwise continuous. 
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Killing X exponentially on the scale of L, we can construct a new Brownian mo- 



tion X onQ. We shall do this using the method of |22 23 1. 

Let K = {Kg, s G M+) denote the right continuous pseudo-inverse of L: 

(1.33) Ks = 'ini{t >0,Lt> s}, s £ M+, 

where — as usual — we make the convention that inf = +oo. The continuity of L 
entails that for every s G M+, Lk^ = s. Clearly, K is increasing, and due to its right 
continuity it is a measurable stochastic process. Fix s G M-|-. It is straightforward to 
check that for every t G M_|_, 

(1.34) {Ks <t} = {Lt > s}. 

Because L is adapted, the set on the right hand side belongs to Tt, and since T is 



right continuous, equality ( |1.34| ) shows that for every s G M+, Ks is a stopping time 
relative to T. We remark that since L only increases when X is at the vertex v, the 
continuity of X implies that for every s G M+ we get X{Ks) = v on {Ks < +oo}. 
On the other hand, we shall argue below that L a.s. increases to +oo, so that we get 
X{Ks) = V a.s. for all s G M+. 

Let /3 > 0. Bring in the additional probability space (IR+, i3(]R+), P/j), where Pg 
is the exponential law with rate /3. Let S denote the associated coordinate random 
variable S{s) = s, s e M+. Define 

We extend X, L, K, and S in the canonical way to these enlarged probability spaces, 
but for simplicity keep the same notation for these quantities. 
Set 

(1.35) C/J = mf{t>0, Lt>5}, 

and observe that since K is measurable we may write C/3 = Ks- Thus as above we 
get X{(is) = V. Define the killed process 



(1.36) Xt 



t < C/3, 

Since this prescription for killing the process X via the (PCHAF) L is slightly dif- 
ferent from the method used in [4,^40J, we cannot directly use the results proved there 
to conclude that the subprocess X of X is still a strong Markov process. However, 
it has been proved in p5] Appendix A] that the strong Markov property is preserved 
under this method of killing, i.e., X is a strong Markov process relative to its natural 
filtration (actually relative to a larger filtration, but we will not use this here). Now 



we may employ the arguments in section 2.7 of 1 22 1, or in section III.2 of 1 35 1 to con- 
clude that X is a strong Markov process with respect to the universal right continuous 
and complete augmentation of its natural filtration. 

It is clear that X has a.s. right continuous paths which admit left limits, and that its 
paths on [0, C/?) are equal to those of X, and thus are continuous on this random time 
interval. Moreover, it is obvious that for every ^ G 5°, we have P^{Ci3 > Hy) = 1. 



BROWNIAN MOTIONS ON STAR GRAPHS 



11 



Therefore, up to its hitting time of the vertex, X is equivalent to a Brownian motion 
on the edge to which ^ belongs, because so is X. Altogether we have proved that — 
under the hypothesis that L is a PCHAF, which will be argued below i n all cases that 
we consider — X is a Brownian motion on Q in the sense of definition |l.l| 

There is a simple, useful relationship between the resolvents R and R of the pro- 
cesses X and X, respectively. Recall our convention that all functions f on Q are 
extended to ^ U {A} by /(A) = 0. 

Lemma 1.11. For all A > 0, / G B{g), ^ G 

(1.37) RxfiO = RxfiO - exiO 4 (e"^^") Rxf{v) 

holds true, where e\ is defined in equation ( |1.15| ). 

Proof. For A > 0, / G B{g), ^ G 5 

Rxf{0 = E^{fj e~^'f{Xt)dt) 

= RxfiO - (e-^«^ e~^'f{Xt+c,) dt) . 
By construction, the last expectation value is equal to 

/3 e-^^y^ e-^'E^(^e-^''^ f{Xt+K.))dtds, 

where we used Fubini's theorem. Consider the expectation value under the integrals, 
and recall that for fixed s G M4-, Kg is an J^-stopping time, while X is strongly 
Markovian relative to T. Hence we can compute as follows 

= E^{e-^^^Ex^K.){f{Xt))) 
= E^{e~^^^)E,{fiXt)), 
where we used the fact that a.s. X{Ks) = v. So far we have established 

RxfiO = RxfiO - E^{e'^<f^) Rxfiv). 

In order to compute the expectation value on the right hand side, we first remark that 
because L is zero until X hits the vertex for the first time, we find that for given 
s G M_|_, Kg > Hy, and therefore Kg = + Kg o 6^^. Hence, and again by the 
strong Markov property, 

= E^[e->^^^)E,{e-^^^), 
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Integrating the last identity against the exponential law in the variable s, we find with 
formula ( |1.15[ ) 

and the proof is finished. □ 



Remark 1.12. Formula (1.37 1 is quite useful, because if the resolvent of X is known, 



then — in view of equation ( 1.15l — it reduces the calculation of Rx to the computa- 



tion of the Laplace transform of the density of C/3 under Py. 

2. The Walsh Process 

The most basic process — which on a single vertex graph plays the same role as a 
reflected Brownian motion on the half line — is the well-known Walsh process, which 
we denote by W = {Wt, t > 0). It corresponds to the case where the parameters a 
and c in the boundary condition ( |1.1[ ) both vanish. This process has been introduced 



by Walsh in 1 38 1 as a generalization of the skew Brownian motion discussed in |20 
Chapter 4.2] to a process in M? which only moves on rays connected to the origin. 

A pathwise construction of the Walsh process in the present context is as follows. 
Consider the paths of the standard Brownian motion B = [Bt, t > 0) on M, and 
its associated reflected Brownian motion \B\ = {\Bt\, t > 0), where | • | denotes 
absolute value. Let Z = {t > 0, Bt = 0}. Then its complement Z'^ is open, and 
hence it is the pairwise disjoint union of a countable family of excursion intervals Ij = 
{tj,tj^i), j e N. Let R = {Rj, j G N) be an independent sequence of identically 
distributed random variables, independent of B, with values in {1, 2, . . . , n} such that 
Rj, j G N, takes the value k € {1, 2, . . . , n} with probability G [0, 1], J2k '^k = 
1. Now define Wt = v if t £ Z, and if t G Ij, and Rj = k set Wt = {k, \Bt\). In 
other words, when starting at ^ G Q°, the process moves as a Brownian motion on the 
edge containing ^ until it hits the vertex at time Hy, and then W performs Brownian 
excursions from the vertex v into the edges Z^, k G {1,2,..., n}, whereby the edge 
Ik is selected with probability Wk- 



As for the standard Brownian motion on M (cf. subsection 1.2 1, we may and will 
assume without loss of generality that W has exclusively continuous paths. 

Walsh has remarked in the epilogue of |38|, cf. also ||2J, that it is not completely 



straightforward to prove that this stochastic process is strongly Markovian. A proof 
of the strong Markov property based on Ito's excursion theory [TS] has been given 
in [36, 37]. A construction of this process via its Feller semigroup can be found in |(2| 
(cf. also the references quoted there for other approaches). 

Next we check that the Walsh process has a generator with boundary condition at 
the vertex given by ( |1.1| ) with a = c = 0. Let / G V{A'^). At the vertex v Dynkin's 
form for the generator reads 



(2.1) A'^f{v) = lim 



Ey(^f{{X{H:^^)))-f{v) 



40 Ey{H^, 
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where H^^ is the hitting time of the complement of the open ball B^{v) of radius 
e > around v. 



Lemma 2.1. For the Walsh process Ey{H^^) = e^. 

Proof. Since by construction W has infinite lifetime, H^^ is the hitting time of the 
set of the n points with local coordinates {k, e), k = I, . . . ,n. Therefore, by the 
independence of the choice of the edge for the values of the excursion, it follows that 
under the stopping time i?^^ has the same law as the hitting time of the point e > 
of a reflected Brownian motion on M+, starting at 0. Thus the statement of the lemma 
follows from equation ( |1.6[ ). □ 

From the construction of W we immediately get 

n 

E 

k=l 

'~ we 



EJf{W{H:^J)]=}_^Wkfk{e), 



with the notation fk{x) = f{k,x), x G M+. Inserting this into equation ( |2?T 
obtain 

n 

A^f{v) = lime-2 J^u;fc(A(e) - f{v)), 
k=i 



and since f'{vk) exists (cf. lemma [L3| ) it is obvious that this entails the condition 

n 

(2.2) Xl^'^^'^^'^) = 0- 



k=l 

For later use we record this result as 

Theorem 2.2. Consider the boundary condition with a = c = 0, and b € 
[0, 1]". Let W be a Walsh process as constructed above with the choice Wk = 6fc, 
k € {1,2,..., n}. Then the generator ofW is 1/2 times the second derivative on 
Q with domain consisting of those f £ Cq{Q) which satisfy condition (1.1b i. 

For the remainder of this section we make the choice a = a = 0, Wk = b^, 
/c e {1,2, inlfTT]). 

Next we compute the resolvent of W. Let A > 0, / E Co{Q), and consider first 
^ = V. Without loss of generality, we may assume that W has been constructed 
pathwise from a standard Brownian motion B as described above, and that B is as in 
subsection 1 .2 Then we get 

n 

Ev{f{Wt)) = Y,bmE^{fm{\Bt\)). 

m=l 

Hence we find for the resolvent of the Walsh process 
(2.3a) Rxfiv)= I r^iv,v)f{v)dr^ 
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with resolvent kernel r^(f , r]), t] ^ Q, given by 

/2A 



^ -V2\d{v,ri) 

(2.3b) r^iv,v) = Yl 



m=l 



and where the integration in (2.3ai is with respect to the Lebesgue measure on Q. 



Now let ^ £ Q. We use the first passage time formula ( |1.26| ) together with formu- 
lae ( |1.15| ) and ( |2.3| ), and obtain 



Lemma 2.3. The resolvent of the Walsh process on Q is given by 

(2.4a) R^fiO = [ r^{C,q)f{r])dv, A > 0, ^ G / G 3(0), 

Jg 

with 

n 

(2.4b) r5:'(C,7?) = r,(e,r?)+ J] 6^,^(6 5,^ -= eA,™(r?), 

A:,m=l 

(2.4c) = 2i(;„ - 5km, 



where r\ is defined in equation ( l.2A\ , and where e\^k, ^x,m denote the restrictions of 



e\ (cf. ( |1.15| )) to the edges Ik, Im respectively. 

Remark 2.4. The matrix S"" = (5^^, k,m = 1,... , n) is the scattering matrix 
as defined in quantum mechanics. We briefly recall its construction in the present 



context, for more details the interested reader is referred to ||28j. is obtained from 
the boundary conditions at the vertex v in the following way. Consider a function 
f on Q which is continuously differentiable in ^° = ^ \ {v}, and such that for all 
k = 1, . . . ,n the limits 

Fk = fivk)= lim /(O 

F'k = f'{vk) = ^ lim /'(O 

exist. Define two column vectors F, F' G C", having the components Fk and F^, k = 
1 , . . . , n, respectively. Furthermore, consider boundary conditions of the following 
form 

(2.5) AF + BF' = 0, 

where A and B are complex n x n matrices. The on-shell scattering matrix at energy 
E > Ois defined as 

(2.6) Sa,b{E) = -{A + iVEB)-\A - iVEB), 

which exists and is unitary, provided the nx2n matrix {A, B) has maximal rank (i.e., 
rank n) and AB'^ is hermitian. These requirements for A and B guarantee that the 
corresponding Laplace operator is self-adjoint on L'^{Q) (with Lebesgue measure). 



Observe that under these conditions the boundary conditions (2.5 1 are equivalent to 



any boundary conditions of the form CAF + CBF' = where C is invertible. Also 
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Sca,cb{E) = Sa,b{E) holds true. For the Walsh process at hand, concrete choices 
for A and B are given by 



A" 



Then ( |2.5| ) is the condition that / is actually continuous at the vertex v, i.e., f{vk) = 
f{vm), k, m = 1, . . . ,n, and that ( |2.2| ) is valid (with Wk = bk, k £ {1,2, . . . , n}). 
Obviously, {A^, B^) has maximal rank. However, A'^{B^)^ is hermitian if and only 
if all bk are equal (i.e., bk = 1/n, k = 1, . . . ,n). Nevertheless, (2^1 exists also in the 
non-hermitian case, and Sj^w (E) = holds for all > due to the relations 
j^wgw ^ _j^w^ Qwgw ^ addition, the following relations are valid: 



/o 





. 


. 


0\ 




Ai 


b2 


63 • 


• bn-l 


bn\ 


1 


-1 


. 


. 













. 











1 


-1 . 


. 













. 














1 . 


. 





, B"" = 






























. 













. 


. 1 


-V 




^0 





. 





0/ 



(2.7) 
(2.8) 



det 5^ 



(-1) 



n+l 



Furthermore, S'^ is a contraction, and the associated Laplace operator is m-dissipative 
on L^{G) since trivially lm{AB'^) = 0, cf. Theorem 2.5 in [26|. When all b^ are equal, 
such that A^(B'^)^ = 0, then is an involutive, orthogonal matrix of the form 



(2.9) 



-1 + 2Pn. 



Pn is the matrix whose entries are equal to 1/n. Pn is a real orthogonal projection. 



that is Pn 



P^ 



P* = p2 jj jjj^Q ^j^jjj dimRanP„ = 1. The 

relation ( 2.4b 1 giving the resolvent in terms of the scattering matrix is actually valid 
in the more general context of arbitrary metric graphs and boundary conditions of the 
form ( |2.5| ), see | |26„29J . 

It is straightforward to compute the inverse Laplace transform of the right hand 
side of formula ( |2.4b| ), and this yields the following result. 

Lemma 2.5. For t > 0, ^, rj £ G the transition density of the Walsh process on Q is 
given by 



(2.10) =p^(t,e,r/) + 

(2.11) 



n 

E 

fc,m=l 



ll^ (^) 2wm g{t,d^{i,7])) li^ (ri) , 



k,m=l 

p{t, ^, rf) is defined in equation ( |1.17 1, p^{t, ^, rj) in equation ( 1.21 1, g is the Gaufi- 
kernel ( |1.4| ), and dy is defined in equation ( |1.16| ). 
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Remark 2.6. Alternatively p'^{t, rj) can also be written as 
(2.12) 



+ E l'^(^) / PdHv&ds)S}:^g{t-s,div,v))hM- 
k,m=i 



Even though this formula appears somewhat more complicated than ( |2.11[ ), it exhibits 



the role of the scattering matrix S"", that is, it describes more clearly what happens 
when the process hits the vertex. 

3. The Elastic Walsh Process 

In this section we consider the boundary conditions ( |1.1| ) with < a < 1 and c = 
0. The corresponding stochastic process, which we will denote by W^, is constructed 
from the Walsh process W of the previous section in a similar way as the elastic 
Brownian motion on M+ is constructed from a reflected Brownian motion (cf., e.g., 
|[T9|, 120} Chapter 2.3], ||23l Chapter 6.2], |[22) Chapter 6.4]). 

In more detail, the construction is as follows. Consider the Walsh process W as 
discussed in the previous section. We may continue to suppose that W has been 
constructed pathwise from a standard Brownian motion B, as it has been described 
there. But then the local time of W at the vertex, denoted by L'^ , is pathwise equal 
to the local time of the Brownian motion at the origin (and we continue to use the 



22 23 35 1) that L 



for the construction 



normalization determined by ( |1.8| )). It is well-known (e.g., |20 
has all properties of a PCHAF as formulated in subsection 1.5 
of a subprocess by killing W at the vertex. We continue to denote the rate of the 
exponential random variable S used there by /? > 0. Let be the subprocess so 



obtained. In particular (cf. 1.5 1, is a Brownian motion on Q, and in analogy with 
the case of a Brownian motion on the real line we call this stochastic process the 
elastic Walsh process. We write q for the lifetime of (i.e., for the random time 



corresponding to in subsection 1.5 1 



We proceed to show that the elastic Walsh process has a generator with 
domain D(A^) which satisfies the boundary conditions as claimed. In other words, we 
claim that there exist a G (0, 1) and 6^ € (0, 1), G {1, 2, . . . , n}, with a + ^j^. bk = 
1, so that for all / G V{G), 

n 

(3.1) af{v) = Y,bknvk) 

k=i 

holds. To this end, we calculate A '^fjv ) in Dynkin's form. We shall use a notation 
similar to the one used in subsection 



1.5 



Namely, let and Ey denote the probability 
and expectation, respectively, on the probability space extended by (M+, i3(M-|-), Pg), 
while the corresponding symbols without " are those for the Walsh process without 
killing. 

For e > and under let ^ denote the hitting time of the complement B^{vY 
of the open ball B^{v) of radius e > with center v. Then ^ = H^^ A C/3,o> 
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where as before H^^ is the hitting time of B^{vY by the Walsh process W. (Note that 
B^{vY contains the cemetery A.) We find 

n 

(3.2) E,(^f{W'{H:^,))) = i^wMe)) < C/3,o). 

k=l 

The probability in the last expression is taken care of by the following lemma. 
Lemma 3.1. For all e, /3 > 0, 

(3-3) n(i?.".<C/3,o) = ^-^. 

Proof. We may consider the Walsh process W as being pathwise constructed from 
a standard Brownian motion B on the real line as in the previous section, and we 
shall use the notations and conventions from there. Then it is clear that under Py 
and under P^, H^^ has the same law as the hitting time of the point e in R+ by the 
reflecting Brownian motion \ B\ under Qq, that is, as H^_^ ^| of the Brownian motion 
B itself under Qq. Let K"^ denote the right continuous pseudo-inverse of L^. For 
fixed s G M+ we get 

{Kf < HI,} = {L^{Hl,) > s}. 

Hence 

P,{Kf < i/- ) = P.{L^{H:^:) > s) 

= Qo(i^(^f_,,,|)>«). 

In appendix B of | (25| is shown with the method in p2{ Section 6.4] that under Qq the 
random variable L^H^_^ ^y) is exponentially distributed with mean e. So we find 

We integrate this relation against the exponential law with rate /3 in the variable s, and 
obtain 

POO 

Pv (C/3,0 > Ke) = 1 - /3 / e-^^ P, {KY < H^,,) ds 

Jo 

1 

~ l + e/3" 

We used the fact that due to the continuity of the paths of W we have C/3,o / ^ 
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We insert formula ( |3.3| ) into equation ( |3.2| ), and obtain 

A''f{v) = lim ^ 



40 E,{H-^, 
lim 



40 E,{H-^, 



1 " 



fiv] 



fiv] 



lim 



40 E,iH^^,)l + eP 



k=l 
n 

k=l 



fk{e)-fiv) 



2. 1 1. Since the last limit and 



Obviously E^{H^J < E^{H^J = (of. lemma 
f'{vk), k G {1,2,. exist and are finite, we get as a necessary condition 



(3.4) 



Y.Wkf'ivk)-f3fiv) = 0. 



k=l 



Thus we have proved the following theorem. 

Theorem 3.2. Consider the boundary condition ( |1.1[ ) with a G (0, 1), h G [0, 1]", 
and c = 0. Set 

h 



(3.5) 



-, k = 1, . . . ,n, P 



— a 1 — a 

and let be the elastic Walsh process as constructed above with these parameters. 
Then the generator of is 1/2 times the second derivative on Q with domain 
consisting of those f G Cq{Q) which satisfy condition ( l.lbl. 

Remark 3.3. Note that condition ( |l.la| ) entails that if Wk and /3 are defined by ( 3.5 1 
then Wk G [0, 1], A; = 1, ... , n, J2k ~ ^' /3 > 0. Therefore the choice ( |3.5[ ) is 
consistent with the conditions on these parameters required by the construction of the 
elastic Walsh process W^. 

Next we compute the resolvent of the elastic Walsh process. As a byproduct 
this will give another proof of theorem 3.2 Moreover, it will provide us with an 
explicit formula for the scattering matrix in this case. In contrast to the calculations 
in pO| Chapter 2.3], p3\ Chapter 6.2] for the classical case with Q = M+, we do 
not use the first passage time formula ( |1.26| ), but instead we use formula ( |1.37[ ). This 
simplifies the computation considerably. 

Let / G Co{Q), A > 0, and ^ G ^. In the present context formula ( |1.37[ ) reads 

where is the resolvent of the Walsh process without killing, and ex is defined 
in ( 1.15 I. The Laplace transform of the density of q under is readily computed: 

Lemma 3.4. For all A, /3 > 0, 



/3- 



2A 
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Proof. As remarked before, we may consider L"^ to be equal to the local time at the 
origin of the Brownian motion B underlying the construction of W, and therefore the 
analogous statement is true for the right continuous pseudo-inverse of L'^ . As 
above let denote the right continuous pseudo-inverse of (cf. O). Then for 

s G M. , 



EJe 



2Xs 



where we used lemma [177] Hence 
which proves the lemma. 



□ 



With lemma 3.4 we obtain the following formula 



(3.6) 



/3 



/3 + V2A 



Note that R^f is in the domain of the generator of the Walsh process, and therefore 
satisfies the boundary condition (2.2 1: 

n 

Y,MRxfy{vk)=o. 

k=l 

On the other hand, we obviously have e'^{vk) 
Thus with J2k=i If fc = 1 we find, 

n 

Y,wu{Rlf)'{vu) = p 



2A 



fc=i 



/3 + V2A 



2A for all k G {l,2,...,n}. 



while equation ( |3.6| ) yields for ^ = 

Rlfiv 



2A 



R'^fiv)- 



/3 + V2A 

The last two equations show that for all / G Co(t/), A > 0, we have 

n 

Y,MRifyM=(3Rif{v). 

k=l 

Since for every A > 0, i?^ maps Co{Q) onto the domain of the generator of W^, we 
have another proof of theorem [3^ 

Upon insertion of the expressions for the resolvent kernels of the Walsh process, 
equations (23 1, and (24i, with the same notation as in lemma 2.3 we immediately 
obtain the following result: 
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Lemma 3.5. For X > 0, i] £ Q the resolvent kernel of the elastic Walsh process 
is given by 



(3.7a) ri(e, r?) = rf v) + Yl ^^''^(^) 



1 



k,m=l 
n 



/3 + V2A 



(3.7b) 



fc,m=l 



w/f/i f/ie scattering matrix 



(3.7c) 5|^(A) = 2 



2A 



Wr. 



X > 0, k,m £ {1,2, .. . ,n}. 



P + V2X 

Remark 3.6. Note that in contrast to the case of the Walsh process, this time the scat- 



tering matrix is not constant with respect to A > 0. Also, when (3 = 0, formula ( |2.4c[ ) 
is recovered, as it should be. In analogy with the discussion in remark [2!4] the bound- 
ary conditions for the elastic Walsh process is given by the matrices 



/o 





. 


. 


M 




( Wl 


W2 


W3 . 


• Wn-1 


Wn\ 


1 


-1 


. 


. 













. 











1 


-1 . 


. 













. 














1 . 


. 













. 













. 


. 1 


-V 









. 





0/ 



such that 



5^(A) = Sa-,b-{E = -2A) 

= -(^^ + \/2AS^)-i(A" 



2XB^). 

Observe that fox k, m G {1, . . . ,n} the matrix element 'S'^„(A) of the scattering 
matrix is obtained from the resolvent kernel as 

5|^(A) = %/2A lim (ri(e,7?)-rA(e,r?)), 

where the limit on the right hand side is taken in such a way that ^, r/ converge to v 
along the edges 1^, Im respectively. S^^{X) in turn fixes the data Wm and /3, e.g., via 
the behavior for large A, that is from the behavior at "large energies" 

Wm = ^(4m+ lim 5'fc„(A')), for all A; G {l,2,...,n}, 

Z A'too 

5fcm + limA'too5'|^(A') 



and 

P = V2X 



(A) 



1 I , for all A, and all A;, m G {1, 2, . . . , n}. 
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Alternatively, the data can be obtained from the small A behavior, that is the threshold 
behavior, of the scattering matrix, since from 

U) 1 

= ^iVi^ {StmW + Skm) for all A: G {1, 2, ... , n}, 

P Zy Za 

we obtain 

1 " 

r ' = ^ (E ^fc-(^) + l) A: G {1, 2, ... , n}, 

m=l 



and therefore 

limA4oA-i/2(5^^(A)+4„) 



lim,^oA-i/2(^^5e,^(A) + l 



for all k, k' £ {1,2, ...,n}. 



Furthermore we remark that in the context of quantum mechanics in the self-adjoint 
case Wk = 1/n, k = 1, . . . ,n, the boundary conditions of the elastic Walsh process 
are interpreted as the presence of a 5-potential of strength (3 at the vertex. 

In order to compute expressions for the transition kernel of the elastic Walsh pro- 
cess, we use the following two inverse Laplace transforms which follow from formu- 
lae (5.3.4) and (5.6.12) in [9] (cf. also appendix C in [25]) (A > 0, t > 0, x > 0): 



(3.8) ^ eoidt)-^{-^^-L^'^/^e.ic{p,l^^))dt, 



Then the inverse Laplace transform of the scattering matrix S**^ is given by the follow- 
ing measures on M+ : 

Skmidt) = C^Wm - 5km) ^o{dt) - 2WmP dt 

(3.10) ^ ^ 

with fc, m G {1, 2, . . . , n}. Moreover, for t > 0, x > 0, let us introduce 

(3.11) g^,(^t,x)=g{t,x)-^ef'^+^''/'evic[-^^+P^-). 

Lemma 3.7. For t > 0, ^, rj £ Q, the transition density ]f of the elastic Walsh process 
is given by 

n 

(3.12) p'{t,i,i^)=p^{t,i,i^)+ h,{i)2wmgp,o{t,d,{i,ii))li^{ii), 

fc, 771=1 
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and alternatively by 

n „i 

(3.13) k,m=l -^0 

X {4m * 9{-, d{v, 7]))^ {t-s)^ l/^(ry), 

where * denotes convolution. 

4. The Walsh Process with a Sticky Vertex 

In this section we construct Brownian motions on Q with a = in the boundary 
condition ( |1.1[ ). 

Consider the Walsh process 1^ on ^ from section [2] together with a right continu- 
ous, complete filtration T"^, relative to which it is strongly Markovian. Furthermore, 
we denote its local time at the vertex v by (cf. section|3]). 



Again we follow closely the recipe given by Ito and McKean in 1 19 1 (cf. also |23 
Section 6.2]) for the case of a Brownian motion on the half line. For 7 > introduce 
a new time scale r by 

(4.1) t + jL"^, t>0. 

Since is non-decreasing, is strictly increasing. Moreover, we haver- 1(0) = 
and limt_s>+oo T^^{t) = +00, which implies that r exists, and is strictly increasing 
from R+ onto M+, too. As is shown in |j23] p. 160], the additivity of entails the 
additivity of r on its own time scale, i.e.: 

Lemma 4.1. For all s, t >0, a.s. the following formula holds true 

(4.2) T{s + t)=T{s) + T{t)o9,^,y 

It is easily checked that for every t > 0, T{t) is an J^"'-stopping time, and since r 
is increasing, we obtain the subfiltration = {Tf, t > 0) of defined by = 
J^^^y t £ R+. Moreover, we set J"^ = cj(J7, t £ M+) and J"^ = a{J^f, t € M+), 
and find C . Standard calculations show that the completeness and the right 
continuity of entail the same properties for T^. (For details of the argument in 
the case where Q = M+ we refer the interested reader to section 3 of |p5|.) 

Define a stochastic process on Q, called Walsh process with sticky vertex, by 

(4.3) Wl = W,^ty teM+. 

Observe that when W is away from the vertex, L"' is constant, and therefore in this 
case grows with rate 1. On the other hand, when W is at the vertex, grows 
faster than with rate 1, and therefore r increases slower than the deterministic time 
scale 1 1— )• t. Thus W'"^ "experiences a slow down in time" until W has left the vertex. 
In this heuristic sense the vertex is "sticky" for W^, because it spends more time there 
than W. 
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Note that because has continuous paths, and therefore also r are pathwise 
continuous. Consequently, W'"^ has continuous sample paths. Since W has contin- 
uous paths, it is a measurable process, and hence for every t > 0, W^(^f^ is -^^(j)- 
measurable, that is, is ^"'^-adapted. Set Of = 9^(ty With the additivity \A2\ of r 
we immediately find 

(4.4) W!oei = W!^t, s, tGM+. 

Thus 9^ = {6f,t€ M+) is a family of shift operators for W^. 

Next we show the strong Markov property of relative to following the 
argument briefly sketched in section 6.2 of [23] for the case Q = ]R_|_. First we prove 
the simple Markov property of with respect to T^. To this end, let s, t > 0, 
^ G ^, and C E 8(0). Then we get with K^ 



= Pw^,,, {W^is) e C) 
= Pws{W!eC), 

where we used the strong Markov property of W with respect to T^. As a next 
step we prove that has the strong Markov property for its hitting time of the 
vertex. By construction, and W have the same paths up to the hitting time of 
the vertex, and in particular Hf is also the hitting time of the vertex by W, that is, 
Hf = Hy. Moreover, since L^{Hy) = 0, we get that T^^{Hy) = = T{Hy), as 
wen as 9^{Hy) = d{Hy). Assume now that t > 0, ^ G a, and C G B{g). Then on 
{Hy < +00} we can compute with the strong Markov property of W as follows 

PdWt+H. e C I -F]^J = PdWf o 9j,^, G C I T^J 
= P.(H^,(t) gC) 

= Py{Wi G C). 

It is readily checked that Tfj C , and therefore we get in particular the strong 
Markov property of with respect to Hf, = in the form 

(4.5) P^ (W,\hs e C I ^I, J = {Wi G C) . 

Finally, with the strong Markov property of the standard one-dimensional Brownian 



motions on every edge and the strong Markov property ( |4.5| ) just proved we can apply 
the arguments similar to those in [20 , Section 3.6] to conclude that is a Feller 
process. Hence it is strongly Markovian relative to the filtration J^*. 

By construction, is up to time equivalent to a standard one-dimensional 
Brownian motion, and it has continuous sample paths. Hence, altogether we have 
shown that W'"^ is a Brownian motion on Q in the sense of definition |1.1| 

Now we want to compute the generator of W^, and first we argue that v is not 
a trap for W^. To this end, we may consider W as constructed from a standard 
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Brownian motion B as described in section [2l Let Z denote the zero set of B. Given 
s > we can choose to ^ s in the complement Z'^ of Z. Consider t = r~^(to), 
i.e., t = to + jL^g- Obviously t > s, and r(t) G Z'^. Therefore -Bt(<) / 0' ^^id 
consequently = 7^ v. 

Theorem 4.2. Consider the boundary condition ( |1.1| ) with a = 0, c G (0, 1), and 

h G [0, 1]". Set 

bk c 

(4.6) Wk = - ,k = l,...,n, 7 = - , 

1 — c 1 — c 

let be the sticky Walsh process as constructed above with these parameters. 
Then the generator of is 1/2 times the second derivative on Q with domain 
consisting of those f £ Cq{Q) which satisfy condition ( |l.lb[ ). 

Before we prove theorem [4!2| we first prepare two preliminary results. Let e > 0, 
and let ^ denote the hitting time of the complement of the open ball Be{v) with 
radius e and center v by W^. Recall that H^^ denotes the corresponding first hitting 
time for the Walsh process W. 

Lemma 4.3. P^—a.s., the formula 

(4-7) i?,% = i7- +7i^|^»^ 

holds true. 

Proof. Let W, and therefore also W^, start in the vertex v. Since and W have 
continuous paths with infinite lifetime we have for all 7 > 

= inf{t > 0, d{v,Writ)) = e}, 

and in particular for 7 = 0, 

H^^, = M{t>0, div,Wt) = €}. 
Moreover, as argued above, both infima are a.s. finite. Set 

Then T{a) = H^^, and therefore 

d{v,W^)=d{v,W^^^)) 
= d{v,WH^J 
= e. 

Consequently we get ^ ^ ^- To derive the converse inequality we remark that 

e = d{v,Wl,,) 
= d{v,W,^Hi^,)), 

which implies 
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Since r is strictly increasing tiiis entails 



and the proof is finished. 
Corollary 4.4. For every 7 > 0, 



□ 



(4.8) 

holds. 



Proof. By construction, the paths of W starting in v hit the complement of Bt:{v) 
exactly when the underlying standard Brownian motion B (cf . section [2]) starting at 
the origin hits one of the points ite on the real line. Thus under Py, L'^{H^^) has 
the same law as L^{H^_^ ^|) under Pq. Lemma 1.8 states that under Pq this random 
variable is exponentially distributed with mean e.Tnen equation (4^) follows directly 
from lemmas 14.31 and|2. 11 □ 



Given these results, we come to the 



Proof of theorem [4!2] Let Wk, k = 1, . . . , n, and 7 be defined as in (4^1, and note 
that due to the condition ( l.la i on bk, k = 1, . . . , n, and c, we have Wk G [0, 1], 
k = 1, . . . ,n, Wk = 1, as well as 7 > 0. Hence we can construct the associated 
sticky Walsh process as above. 

Let A'^ denote the generator of with domain 'D{A^). Then we have for / G 
V{A^), A^f{v) = l/2f"{v) (cf. theorem 1.5 I. On the other hand, we can compute 
via Dynkin's formula as follows 



lim 

40 



lim 

£4-0 



E,(^f{W^{H^J))-f{v) 



EJm 



Y.k'^kfkje) - f{v) 
+ 76 



where we used corollary 4.4 Since the directional derivatives of / at t> 

/(e) - f{v) 



f'iVk) 



lim 



diC,v) 



k G {l,2,...,n}, 



exist (cf. lemma [T3] ), we obviously get the boundary condition 

1 



(4.9) 



1 

z 7 ^ — ' 



k=l 



as a necessary condition. Finally, inserting of the values ( |4.6| ) of the parameters Wk, 
k = 1, . . . , n, and 7 into equation (|4.9|) we complete the proof of theorem [42 



□ 



Next we shall compute the resolvent of the Walsh process with sticky vertex. 



Similarly to the alternative proof of theorem 3.2 for the elastic Walsh process, as a 



byproduct we obtain an alternative proof of theorem 4.2 We begin with the following 
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Lemma 4.5. Let X > 0, f e Co{g). Then 

(4.10) l{Rlf)"{v) = ^ (2A(e^,/) - V^fiv)) 



2' ^/2A + 7A 

holds, where 

(4.11) e^(0 =«^fceA(0, i£lk,k = l,...,n, 

and ex is defined in equation ( 1.15| ). 

Proof. Let A"* be the generator of VF"* on Co (^). From the identity A^'i?^ = Ai?^-id, 
and the definition of r we get 



I {RVy'iv) = \E,[j^ e-''{fiWn - f{v)) dt) 

= XE, e-^(^+'^^=^) { f{Ws) - fiv)) {ds + ^dLf)) 
= XE, ( r e-^(^+^^«"') {fiWs) - f{v)) ds) . 



In the last equality we used the fact that only grows when W is at the vertex v. 
By construction of the Walsh process W we have 



EJe-^^'^^f{Ws)-fiv))) 

n 

--^WkEo[e-'^''^fk{\Bs\) - fk{0)] 

k=l 

^ f'OO f'OO I 

= 2^^,/ / e-^^^(/fe(x)-A(0))4±|_e-(-+^)V2^dxdy, 
^ Jo Jo VItts-^ 



where we used lemma 1 .6 We insert the last expression above, and use formula ( |1.31 
This gives 

1 ^ -| poo 

{Rlf)"{v) = 2Xj2^,-= / e-^'^{Ux)-MO))dx 

^ ^ \/2A + 7A Jo 

^ (2A(e^,/)-^/2A/(t;)). □ 



^2A + 7A 

From the identity A^R\ = A i?^ — id and some simple algebra we get the 
Corollary 4.6. Let A > 0, and f G Co(^). T/jew 

(4.12) = -jJ— (2(e^,/) +7/W) 

V 2A + 7A 

holds. 
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Since formula ( |1.27| ) in corollary 1. 10 is valid for the resolvent of every Brownian 
motion on Q, we may use that formula for sum it against the weights Wk, k = 
1, . . . , n, and insert the right hand side of equation ( |4.12[ ). This results in 

1 



k=l 



2A + 7A 



(2A(e^,/)-V2A/(«)), 



and a comparison with formula ( 4.10| ) shows that equation (4^1 holds true for / re- 
placed by i?^/ for arbitrary / G Co{Q). As promised we thus have another proof of 
theorem [4.2[ 

With the help of the first passage time formula we can now provide explicit expres- 
sions for the resolvent R^, its kernel and the transition kernel of W^. Inserting 
the right hand side of equation ( |4.12 1 into the first passage time formula ( 1.26 l, we 
immediately obtain for / € Co(^), A > 0, 

(4.13) Rlf(C) = R^f{0+ ^ eA(0(2(e^,/) + 7/W), UG, 



^2A + 7A 

where R^ is the Dirichlet resolvent ( 1.22| ). Using formula ( 1.23 1 for the kernel of R^ 
together with ( 1.26 ), and ( 1.27 1, we get the following result. 

Corollary 4.7. For ^, rj £ Q, \ > 0, the resolvent kernel r|, of the Walsh process 
with sticky vertex is given by 



(4.14) 



1 



k,m=l 



+ 



2A + 7A 

7 



2A + 7A 



with defined in (1.23 1, and denotes the Dirac measure in v. Alternatively, is 
given by 

n ^ 

riitdr]) = rx{C,v)dv+ Yl Si^{\) ^= ex,m{ri) dr] 



(4.15a) 



k,m=l 



+ 



7 



2A + 7A 



eA(6 ^v{dr]), 



where r\ is defined in equation ( 1.24 1, and 



(4.15b) 



SkmW 



2A 



2A + 7A 



^km- 



Remark 4.8. When all Wm, m = 1, . . . , n, are equal to 1/n, the matrix S^{X) takes 
the form 

S%\) = -1+ ^ P„ 



2A + 7A 



which reduces to ( |2.9| ) when 7 = 0. 5* (A) is unitary for all A < 0. Also the S^{X) 
for different A all commute. As a consequence 5^ (A) has the interpretation of a 
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quantum scattering matrix in the sense of [28 1. More precisely, <S'^(A) stems from 
the Schrodinger operator — A'^, where A'' is a self-adjoint Laplace operator on L'^{Q) 
with boundary conditions of the form (|2.5|) with the choice 



^ = -2(^^(^o) + l), 
(4.16) 1 

2? = -— ==(5^(Ao)-l), 



for any Aq for which \/2Ao + 7A0 7^ 0. We emphasize that the Schrodinger operator 
—A* and the generator ^4'^ of the Walsh process are quite different: Not only do 
they act on different Banach spaces, but also the functions in the intersection of their 
domains satisfy different boundary conditions at the vertex v. As matter of fact, the 
integral kernel of the resolvent (— A^ + 2A)^^ of the Schrodinger operator —A* is 
given by, see Lemma 4.2 in p9}. 



1 1 



k,m=l 

that is — up to a factor 2 — by the right hand side of ( 4.15a| ) without the last term. 
In more detail and with the definition ( |2.6[ ) 

Sa,b{E = -2X) = S%X) 

holds for all A > 0. As a function of k (k'^ = E), is meromorphic in the complex 
/c-plane with a pole on the positive imaginary axis at k^ = 2i/ 7. This corresponds 
to a negative eigenvalue = —4/7^ of — A'^. The corresponding (normalized) 
eigenfunction — physically speaking a bound state — is given as 

1 Ff _-2d(t,,0/7 



2 \ n 

So quantum mechanically the vertex v acts like an attractive potential. We view this 
as a quantum analogue of the stickiness of the vertex v. 

This analogy can be elaborated a bit further by inspecting the associated quantum 
mechanical time delay matrix (see, e.g., p^[5|[2T][3T} 



which in the present context gives 

^^^^ = fe(4+P7) 
So T{k) has zero as an (n — l)-fold eigenvalue plus the non-degenerate eigenvalue 

-27 
A;(4 + fcV)' 

which for 7 > is the signal for a strict quantum delay. Observe that for k — )• +00, 
that is for large energies, the time delay experienced by the quantum particle tends to 
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zero, while for A; — 0, i.e., for low energies, the delay becomes arbitrarily large. From 
the physical point of view, both effects are clearly to be expected. For comparison 
and in contrast to the present stochastic context, in quantum mechanics 7 < is also 
allowed for a meaningful Schrodinger operator and an associated scattering matrix. 



Define for x > 0, 7, t > 0, 



go,j{t,x) = - exp 
7 



(4.17) 

It is not hard to check that 
(4.18) 



2x 2t 

7 T 



eric I —;= H 

7 



2t 



limgo,^{t,x) = g{t,x) 
74,0 



/2Trt 



Moreover, from ||9j eq. (5.6.16)] (cf. also appendix C in |25 1) the Laplace transform 
is 

1 



(4.19) 



2Xx 



2A + 7A 



X > 0. 



Observe that in agreement with ( |4.18[ ) 

Cg{-,x)iX) 



2Xx 



2X 



holds. Now we can readily compute the inverse Laplace transform of formulae ( |4.14| ), 
( 4.15 1, and obtain the following result. 

Corollary 4.9. For t > 0, ^, r] £ Q the transition kernel of the Walsh process with 
sticky vertex is given by 



(4.20) 



m 50,7 ; dy ('^; 1/m (^) drj 

k,Tn=l 



+ 79o,j{t,d{^,v)) ey{dr]) 

where is defined in equation ( |1.21[ ), or alternatively by 
p\t,i,dr]) = p{t,^,ri)dr] 

n 

m go,^{t,dv{S,,rj)^ drj 

(4.21) k,m=l 

- 5kmgit,dv{(,,r])) di]^ 

+ 79o,y{t,d{^,v)) ey{d7]), 
andp{t, a, b) is given in formula ( 1.17[ ). 
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We close this section with some remarks concerning the local time of at the 
vertex v, which also serve to prepare the construction of the most general Brownian 
motion on the single vertex graph Q in the next section. 

Let us define 

(4.22) Ll = KitY i>0, 

where — as before — denotes the local time of the Walsh process at the vertex, 
having (cf. section [3]l the same normalization as the local time of a standard one- 
dimensional Brownian motion (cf. ( |1.8| )). By construction, is pathwise continuous 
and non-decreasing. It is adapted to F^, and a straightforward calculation based on 
the additivity of and formula ( |4.2| ) shows the (pathwise) additivity property 

(4.23) Ll^t = LI + LI oei s,t>0. 

Thus L^ is a PCHAF of {W^ ,F^). Furthermore, t > is a point of increase for 
if only if r(t) is a point of increase for L"', which only is the case if W^-i^f^ is at the 
vertex, i.e., if is at the vertex. Thus, it follows that L^ is a local time at the vertex 
for W^. In order to completely identify it, it remains to compute its normalization, 
and it is not very hard to compute its a-potential (the interested reader can find the 
details for the case Q = M+ in [25J): 

(4.24) EJre-"''dLt)=^=J^ ^-V2ad{^,v)^ a>0,^eg. 



'2a + 7a 

5. The General Brownian Motion on a Single Vertex Graph 

Finally, in this subsection we construct a Brownian motion by killing the Walsh 
process with sticky vertex of section [4] in a similar way as in the construction of the 
elastic Walsh process (cf. section [3]l. realizes the boundary condition \\.\\ in its 
most general form. 

Consider the sticky Walsh process with stickiness parameter 7 > 0, right 
continuous and complete filtration , and local time L^ at the vertex. We argued 
in section |4] that L' is a PCHAF for and therefore we can apply the 

method of killing described in section |1.5[ We bring in the additional probability 
space (M+, ;S(M+), P^) where is the exponential law of rate /3 > 0, and the 
canonical coordinate random variable S. Then we take the family of product spaces 
[Vl^A, {P^, e G G)) of {n,A, {P^, ( e g)) and (M+, e(M+), P^). Define the ran- 
dom time 

(5.1) C/3,7 = mf{t > 0, > 5}. 

Then by the arguments given in section |1.5[ the stochastic process defined by 
= for t £ [0, C/3.7)> and Wf = A for t > (is^j, is again a Brownian motion 
on Q in the sense of definition |1.1| 

Denote by the right continuous pseudo-inverse of L*. Since L"^ is continuous 
(cf. equation < \4.22 1), we get L|^s = r for all r G M^. Recall that the right continuous 



pseudo-inverse of the local time L^ of the Walsh process was denoted by K"^. Then 
we have the following 
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Lemma 5.1. For all 7 > 0, the following relation holds true: 
(5.2) K'^=K';^ +-ir, rEM+. 

Proof. For 7, r € M+ define the random subset 

J^(r) = {t > 0, > r} 

of M+. Since is pathwise increasing, this set is a random interval with endpoints 
K?. and +00. The relation Ll-s = r implies that 

r 

j^(r) = (^^^^+oo). 

In particular, we have Jo = (-^^j +00). Now 

t e J^{r) 4^ Lf = L;f(^) > r 4^ T{t) € Jo(r). 

In other words, J^{r) = [Jo{r)), and therefore = t~^{KI") holds. From 
the definition of r^^ (see equation ( |4.1| )), and the relation = r we obtain for- 



mula d53 □ 



In the proof of lemma 3.4 the Laplace transform of the density of K'^, r > 0, under 



has been determined as A 1— )• exp(— v 2Ar). Hence we have 



P^{Kl" G dl) = e"" I > 0. 

As a consequence we find the 

Corollary 5.2. For r > 0, has the density 

(5.3) P^iKl G dl) = , g-rV2(«-7r) ^ > 

Furthermore, the Laplace transform of the density of under P^ is given by 

(5.4) S,(e-^^^^) =e-(^+^^>, A > 0. 



Remark 5.3. One can use lemma C.l in |25| to check that the right hand side of 
equation (53 1 is indeed the inverse Laplace transform of the right hand side of for- 
mula ( [54l ). 

Observe that = Kg and C/3,0 = Kg. Thus we obtain the 
Corollary 5.4. For all P > 0, ^ > 0, the following equation holds true 
(5-5) C/3,7 = C/3,0 + iS. 

As before, denotes the expectation with respect to Pg, ^ G Q. 
Corollary 5.5. For a// /3 > 0, 7 > 0, A > 0, the following formula holds true 
(5.6a) ^.(e-^^"--) =/3/>(A), 

(5.6b) p{X) ^ 



/3 + V2A + 7A 
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Proof. With corollary |5 .21 and C/3,7 



Kg we obtain 



/3 + V2A + 7A 



□ 



Denote by W the resolvent of . With lemma [TTTT] we immediately find the 

Corollary 5.6. For all f G Co{Q), X > 0, £ G the following formula holds true: 

(5.7) RlfiO = RlfiO - I3p{\) exiO Rlfiv). 

Now it is easy to verify that for appropriately chosen parameters /3, 7, w^, k = 
1, . . . , n, the Brownian motion Wg realizes the boundary condition ( l.lb i. 

Theorem 5.7. Consider the boundary condition ( |1.1[ ), and assume that b is not the 
null vector Set r = a + c G (0, 1), and 



(5.8) 



bk 



-, k = l,...,n, j3 



7 



Let be the Brownian motion as constructed above with these parameters. Then 
the generator ofW^ is 1/2 times the Laplace operator on Q with domain T)[A^) 
consisting of those f € Cq{Q) which satisfy condition (1.1b). 



Proof. As in the previous cases it is readily seen that the definition ( |5.8| ) of the pa- 
rameters 7, (3, Wk, k = 1, . . . , n, is consistent with the conditions used in the above 
construction of W^. 

Let be the generator of with domain D(A^). Since is a Brownian 
motion on Q in the sense of definition |1.1[ it follows from theorem 1.5 that 'D{A3) C 
C^iQ), and that for all / G V{A9), A9f{0 = ^I'^f'iO, C e Let /i G Co(a), 

follows 



A > 0. Then R\h G 'D{A9), and therefore we may compute with equation (|5.7|) as 



7 



^ {Rih)"{v) = 1 {Rlh)"{v)-Pp{\)2\{Rlh){v) 

n 

= ^wk {Rih)'{vk) -Pp{X)jX{Rih){v), 

k=l 

where we used the fact that, since Rj^h is in the domain of the generator A^ of W^, it 
satisfies the boundary condition ( |4.9| ). We rewrite this equation in the following way: 

n 

I {Rlh)"{v) = ^wk {Rihy{vk) + (3V2Xp{X){Rih){v) 
(5.9) ^ k=i 

-/3p(A)(^/2A + 7A) {Rlh){v). 
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Now we differentiate equation ( |5.7[ ) at ^ G Ij^, k = 1, . . . , n, let ^ tend to v along any 
edge Ij^, and sum the resulting equation against the weights Wk, k = I, . . . ,n. Then 
we get the following formula 

n n 

(5.10) Y.'^k {Rih)'{vu) = Y.'"k {Rih)'{vk)+pV2\p{\){Rih){v), 

k=l k=l 

where we used Wk = 1. On the other hand, for ^ = v, equation ( |5.7| ) gives 

(5.11) [Rlh] {v)=p{X){^/2\ + 7A) {Rih) (v) . 

A comparison of equations ( |5.10| ), ( |5.11| ) with ( |5.9| ) shows that we have proved the 
following formula 



(5.12) 



I {Rih)"{v) = Y,w, {Rih)'{v,) - P{Rih){v). 

k=l 



With the values ( |5.8[ ) for /3, 7, and w^, k = 1, . . . , n, it is obvious that / = R^^h 
satisfies equation ( l.lb| ). Since ii^ is surjective from Co{Q) onto the domain of the 
generator of W^, the proof of the theorem is finished. □ 



Let A > 0, / G Cq{G)- Insertion of the right hand side of formula ( 4.13 1 for i?^ 
into equation ( |5.7| ) gives us after some simple algebra the following expression for 

(5. 13) RifiO = R^fiO + P(A) exiC) (2(6^, /) + 7f{v)) , ^gG, 

where R^ is the Dirichlet resolvent, ex is defined in equation ( |1.15| ), in equa- 
tion ( |4.11[ ), and p{X) is as in formula ( |5.6b[ ). From equation ( |5.13 1 we can read off 
the following result: 

Corollary 5.8. For S,, rj ^ G, X > 0, the resolvent kernel of the general Brownian 
motion on G is given by 

n 

rlitdv) = rx{^,v)dr] + ^ ex,k{0'^Wjn p{X} ex,m{ri) dr] 



(5.14) 



fc,m=l 



with as in formula ( 1.23| ), and p is defined in equation ( |5.6b[ ). Alternatively, 
can be written in the following form 

n ^ 

r{itdr]) = rxi(„7])dr]+ V ex,k{C) Sl^{X) ex,m{r]) dv 
(5- 15a) 1^^^^^ V2A 

+ 7p(A)eA(0 (^v{dv), 

where rx is defined in equation <\1.24 1, and 
(5.15b) SlJX) = 2V2Xp{X) 



34 



V. KOSTRYKIN, J. POTTHOFF, AND R. SCHRADER 



In order to invert the Laplace transforms in equations ( |5.14[ ), ( |5.15| ), we define for 
/?, 7 > 0, tlie following function g^^^^ on (0, +oo) x M_|_: 



with G (0,+oo) X M+. The heat kernel g^^^ is discussed in more detail in 

appendix C of | [25| . In particular, it is outlined there that the limits of gp^.y as /3 | 
0, and 7^0, yield the kernels gp^ (equation ( |3.11| )) and 510,7 (equation ( |4.17| )), 
respectively. Moreover, it is proved there that the Laplace transform of gj3^^{-,x), 
a; > 0, is given by 



(5.17) p(A)e-^2Ax^ ^>Q^ 
where p is defined in ( 5.6b| ). Hence we get the 

Corollary 5.9. For i] £ Q, t > 0, the transition kernel of the general Brownian 
motion on Q is given by 

n 

(5.18) + '^h{02w„,gp,^{t,d^{C,ri))li^{r])dri 

k,m=l 

+ 19l3,'y{t,d{^,v)) ev{dr]), 
which alternatively can be written as 



n 



(5.19) k,m=l 

- 5km9{t,dy{S,,r]))^ lu(^) dv 
+ 1 9Pn{^^d{i,v)) ey{drj). 
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